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ABSTRACT 

It  is  shown  here  that  the  permanent  of  an  n x n doubly  stochastic  matrix 
is  not  less  than  e n. 
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SIGNIFICANCE  AND  EXPLANATION 


Let  A be  a square  matrix.  The  permanent  p(A)  of  A is  basically  the 
determinant  of  A where  all  the  summands  appear  with  + signs.  T'.-i  notion 
p (A)  arises  naturally  in  many  combinatorial  settings  where  a count  of  the 
number  of  systems  of  distinct  representatives  of  some  configuration  is  required. 

In  this  paper  we  establish  lower  estimates  of  the  right  order  for  the  permanents  of 
doubly  stochastic  matrices.  Recall  that  A is  called  a doubly  stochastic  matrix 
if  all  the  entries  of  A are  nonnegative  and  each  row  and  column  sum  of  A is 
equal  to  1.  Doubly  stochastic  matrices  appear  frequently  in  probability  and 
combinatorics.  The  result  established  in  this  paper  will  have  various  applica- 
tions in  combinatorics  and  probability,  in  particular  to  Latin  squares  and 


block  designs. 
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I.  Introduction. 

Let  A be  an  n * n matrix  (a^)".  The  permanent  of  A is  defined  by 

U.l>  p(A)  - l n a , 

oeS  i-1  ' ’ 

n 

where  S is  the  symmetric  group  of  order  n.  Let  ft  be  the  set  of  all  doubly  stochastic 
n n 

matrices*  that  is*  the  set  of  all  n * n matrices  A satisfying 

n n 

(1.2)  i °-  I 'ij  * I “ij  - l>  • 

13  j-1  13  i-1  13 

It  is  conjectured  that  any  doubly  stochastic  matrix  satisfies  the  inequality 

(1.3)  p(A)  >_  nl/nn  ( 

with  the  equality  holding  only  for  the  matrix  J all  of  whose  entries  are  1/n.  The  problem 

n 

of  finding  the  minimum  of  p(A)  on  the  set  fi  goes  back  to  van  der  Waerden  [6] . In  fact 

n 

the  inequality  (1.3)  is  comnonly  referred  as  the  van  der  Waerden  conjecture.  This  conjecture 
is  known  to  have  applications  to  certain  combinatorial  problems.  In  this  paper  we  establish 
the  inequality 

(1.4)  p(A)  > e'n,  A £ 0 

— n 

Recall  that  by  Stirling's  formula  nl/nn  */2nn  e n.  The  previous  lower  bound  known 
before  was  1/nt  . See  (3)  for  the  proof  of  this  result  and  for  the  survey  of  the  main  re- 
sults achieved  in  connection  with  the  van  der  Waerden  conjecture.  Our  starting  point  is  the 
inequality  due  to  T.  Bang  (1) 
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2.  Preliminary  results. 


Let  a e S . Define 
n 


(2.1) 


Po  “ (6a(i)j,l 


By  t we  denote  the  set  of  all  matrices  a with  integer  coefficients  such  that  a/m  is 
n,m 

an  n x n doubly  stochastic  matrix: 

(2.2) 


» _ ■ {a|a  - (a.  .)",  a . t SZ  , o/m  t (i  } . 

n,m  1 ij  1 ij  n 


From  the  classical  Birkhoff  theorem  (e.g.  [4])  it  follows  that 

m 


(2.3) 


a - l P , at*  . 

W °k 


Following  T.  Bang  [1]  we  bring  a formula  for  the  permanent  of  A 9 J 


Theorem  2.1.  Let  A be  an  n * n matrix.  Then 


(2.4) 


p(A  • J ) 


. ,.2n  n a.. 

i»U V n - ij 


m**"  at*  i,  j“l 
n pin 


l n a..  Vaijl  . 


Proof.  Consider  a term  in  p(A  • J ) . It  is  of  the  form 

m 


(2.5) 


n a 

i.j-1 


ij 

ij 


Clearly,  each  a^  is  a non-negative  integer.  Consider  the  rows  i,  m+i , . . . , (n-l)m«-i . These 


n a.  , 
-n  n . ij 


rows  contribute  to  the  product  (2.5)  the  term  m n a,  * . Therefore,  J a « m.  In 

j-1  13  j-1  13 

n 

the  same  way  one  shows  that  7 a « m.  Thus,  («,.,>?«•  • Vice  versa,  s<q>pose  that 

ij  Ij  1 “»■ 

at*  . Then  (2.3)  holds.  Let  us  view  A • J as  a block  Mtrix 
n,m  m 


(2.6) 


-3- 


1 


In  the  block  A take  a product  m Ii  a.  , ...  By  multiplying  all  these  terms  to- 
rr  .,10(1) 

i=l  r 

gether  one  obtains  the  expression 


m n 

(2.7)  m"nm  n n a.  . 

...  l a (l) 
r=l  i=l  r 

In  view  of  (2.3)  the  expression  (2.7)  is  equal  to  (2.5).  We  now  compute  the  coefficient  of 
the  term  (2.5).  That  is  we  are  looking  for  the  number  of  different  ways  to  pick  up  a ele- 
ments a /m,...,a  elements  a /m  from  the  matrix  A ® J such  that  any  two  elements  are 

11  nn  nn  m 

not  on  the  same  row  or  column.  We  call  such  a choice  an  admissable  choice.  Let  us  label 
a.,  elements  a../m  by  a f V /m, . . . , a.  . ^ /m.  Assume  that  for  a.  . > 1 the  element  a(k'/m 

i]  1 3 13  13  13  ~ 13 

(1  £ k <_  a„)  sits  in  the  row  i + u(k,i,j)n  and  the  column  j + v(k,i,j)n  in  the  matrix 
A ® J . We  assume  that  no  two  elements  are  placed  in  the  same  position.  Let  us  call  such  a 

choice  of  positions  a configuration.  Two  configurations  considered  to  be  equal  if  for  each 

(k) 

a. . > 1 and  1 < k < a. . the  Dositions  of  the  elements  a. . /m  coincide.  Given  a con- 
13  “ - - 13  - i] 

figuration  one  obtains  a distinct  conf irguration  by  interchanging  rows  (columns)  i and  j, 
where  i : j (mod  n) . It  is  easy  to  see  that  one  can  obtain  any  admissable  configuration  from 
a given  configuration  by  interchanging  appropriate  rows  and  columns.  Obviously,  the  rows 
(columns)  i,  i+n,...,  i + (m-l)n  (1  <_  i £ n)  can  be  interchanged  in  m!  ways.  Thus  accord- 
ing to  what  we  proved  one  has  (m!)2n  distinct  configurations.  Let  us  go  back  to  the  problem 
of  determining  the  number  of  different  ways  to  pick  up  the  elements  a_/m  f°r 

i,j  = l,...,m  (an  admissable  choice).  Clearly  any  configuration  qives  rise  to  an  admissable 

n 

choice.  We  claim  that  to  any  admissable  choice  correspond  n a. .!  distinct  configuration. 

i.3=l  13 

Indeed,  for  this  choice,  we  have  a..(>l)  places  occupied  by  a. ,/m.  In  these  a. . places 

(1)  *ail'  13  13  13 

we  put  a^  /■,..., a^  /m.  This  can  be  done  in  a„t  ways.  Thus,  to  the  given  admissable 

n 

choice  correspond  n a. ,1  distinct  configurations.  Obviously,  to  two  distinct  admissable 
i , j=l  3 

choices  correspond  distinct  configurations.  This  demonstrates  that  the  number  of  diff<  rent 

ways  to  pick  up  a.,  elements  a ,/m,  i,j  = l,...,m,  from  the  matrix  A R J , such  that 
13  13  _ m 

o n 

any  two  elements  are  not  on  the  same  row  or  column,  is  equal  to  (m!)'n(  n a ! ) "*  . The 

. . , lj 

10*1 


proof  of  the  theorem  is  completed. 

The  permanent  of  a non-negative  matrix  can  be  estimated  in  terms  of  the  permanent  of 
A ® J [1] . 


Theorem  2.2.  Let  s be  a positive  integer  and  m = 2s.  Then  for  any  non-negative  square 
matrix  the  following  inequality  holds 

. 1/m 


(2.8) 


p(A)  > [p(A  ® J ) ) 
m 


Proof.  We  prove  first  the  inequality  (2.8)  for  m » 2.  From  (2.4)  it  follows 

n a. 


(2.9) 


p(A  ® J ) - [ n a ij/a  I 

ae*  , i,j-l  13  13 

n,z 


On  the  other  hand 

(2.10) 


2 n n 

p (A)  » l n a n a.... 

o.nts  i-1  10(i)  j-1  *n(j) 

n 


Let  a 6 ♦ According  to  (2.3)  a = p + p * Thus 

n*2  a n 


(2.11) 


ij 


n n n a 

It  s . / j » n ii»  * n s . i • 

i-1  i0(i)  j-1  Jr,<J)  i,j-l  ij 


n aii  2 

Therefore,  the  coefficient  of  the  term  11  a,  . J in  p(A)  is  a positive  integer.  The 

i.j-1  3 

coefficient  of  this  term  in  p(A  ® J^)  never  exceeds  1.  This  establishes  the  inequality 
(2.8)  for  m - 2.  Now  the  general  case  easily  follows  by  induction.  Indeed,  suppose  that 
(2.8)  holds.  Then 


(2.12) 


p(  (A  ® J > X J,)1/2“  < p (A  e JJ1/n  < p(A) 

n a ui 


It  is  left  to  note  that  the  tensor  product  is  associative  and 


(2.13) 

End  of  proof. 


J»  ® J2  ’ 3 2m 


Let  A be  an  n * n doubly  stochastic  matrix.  Suppose  that  the  van  der  Waerden 


•5' 


conjecture  holds.  Then 


[ (nm) !/(nm)  ] 


In  what  follows  we  estimate 


be  a non-negative  matrix.  Then 


we  get  that  the  coefficient  of  the  term  fl  a 


Choosing  8 


Expanding  the  expression 


appears  in  (2.17)  with  the  coefficient 


we  see  that  the  term 


-iSU tl  (a  .!)  . As  the  expansion  of  the  term  (2.17)  contains  only  non-negative  terms 

mn  . , , lj 
m 1,3=1 

from  the  identity  (2.4)  we  deduce  the  inequality  (2.15). 


Thus,  if  A is  a stochastic  matrix  the  inequality  (2.15)  implies 


Note  that  if  A = P then  the  equality  sign  holds  in  (2.18).  Recall  the  well  known  Stirling 


formula  (e.g.  (2,  p . 52 1 ) 


Thus,  if  the  van  der  Waerden  conjecture  holds  by  combining  (2.14)  and  (2.18)  with  the 


Stirling  formula  we  obtain 


(2.20) 


-n 

e 


for  any  doubly  stochastic  matrix  A.  We  shall  prove  the  equality  above  in  the  next  section, 
without  assuming  the  validity  of  the  van  der  Waerden  conjecture.  Of  course,  once  (2.20)  is 
established  in  view  of  Theorem  2.2  one  deduces  the  inequality 


(2.21) 


p(A)  > e ",  A e f!n  . 


Vice  versa,  (2.21)  implies 

(2.22) 


p(A  ® J ) 
m 


l/m 


Combining  (2.22)  with  (2.18)  we  obtain  again  the  equality  (2.20).  This  manifests  the  equi- 
valence of  the  relations  (2.20)  and  (2.21). 


(3.4) 


(nk^)  1 . . . (nk^)  1 


nk  +1/2  nk  +1/2 

es(nk  ) ...  (nk  ) S 

1 s 


iv  n"  iv  Mn  “ n(k  +1/2)  n(k  +1/2) 

(kl!)  <2H)  1 S/2(k.)  1 ...(k>  S 

1 s 


s mn+s/2 
e n 


s mn+s/2 
e n 


(<2II)sk  ...k  ]<n_1)/2  [(2H)s(m-s+l)]  (n_1)/2 

l s 


Here  we  used  the  inequality 


k,  — k > (m-s+1) 
1 s — 


since  kj,...,k  are  positive  integers  which  sum  up  to  m.  Clearly  the  relations 


(3.5) 


l+s/2nm 


lim 


= lim 


s/nm  l+s/2nm 
e n 


(e"(2nm/s)  <n_1)/2]s/nn’  [ (2!I)  s (m-s+1)  ] (n_1)/2nln 


prove  the  lemma. 

Lemma  3.2.  Let  p^,...,pr  be  non-negative  integers  such  that 


(3.6) 


Pj+.-.+p^  = nm  . 


Then  there  exist  non-negative  integers  ql(...qr  with  the  following  properties 


(3.7) 


ql+‘”+qr  = 1,1 


(3.8) 

and 


I Pj  ~ n9jl  < n,  j = 1, . . . ,r 


(3.9) 


(nm)  I < / n \r  1 (nm)  1 

ll...pr!  - I O (nqi)!...(nqr)! 


Proof.  By  rearranging  the  indices  we  may  assure  that 

Pj  $ o (mod  n) , j = l,...,k  , 


(3.10) 


Pj  i 0 (mod  n) , j ” k+l,...,r 


In  case  that  k » 0 (3.9)  is  trivial.  Assume  that  k >_  2.  Without  loss  of  generality  we 


: _ 


-9- 


assume  that 


Then 


Noting  that  the  function  (nx+a)/(nx+b)  is  increasing  on  (0,“)  if  0 £ a <_  b we  obtain 


Recalling  the  well  known  inequality 


one  deduces 


In  case  that  t < t 


given  by  (3.17)  by  d -t_  and  n(q.+l) 


respectively  and  (3.18)  is  still  valid.  Note  that  p*+...+p 


are  not  divisible  by  n.  Continuing  the  procedure  above  we  obtain  the  in' 


equality  (3.9) 


Proof.  Recall  that 


According  to  Lemna  3.2  there  exist  positive  integers  q 


such  that  (3.7)  - (3.9) 


As  o.  < 1 we  also  have 


For  a fixed  r integers  q. 


q we  can  have  at  most  (2n-l)  types  of  r integers 


p^ satisfying  (3.8).  Thus,  by  using  the  inequality  (3.24)  for  each  summand  appearing 

in  (3.21)  we  establish  the  leirma. 

We  are  now  ready  to  prove  our  main  result. 

Theorem  3.1.  Let  A be  an  n x n doubly  stochastic  matrix.  Then 

(3.25)  lim  p ( A ® J ) = e n . 

m«» 

Proof.  From  the  classical  Birkhoff  theorem  (e.g.  (41)  it  follows  that 

r r 

(3.26)  A = l up  , a > 0,  j = l,....r,  £ 0=1  . 

j=l  3 °j  3 j=l  3 

Without  loss  of  generality  one  may  assume  that  {a  is  a decreasing  sequence,  i.e.  (3.19) 


holds.  We  claim  that 


p(A  ® J ) > 


nk^  nkf 
“l  ••■“rr 


m kj+...+k  =m  (k^!)  . . . ( k^. ! ) 


Indeed,  let  k^,...,kr  be  r non-negative  integers  which  sum  up  to  m.  Consider 

i3-m’  ' ■ L ""S  • ■ 

From  the  expansion  (2.4)  for  the  permanent  of  p(A  ® Jm>  it  follows  that  one  has  a term  of 
the  form 

,_.,2n  r 8ij 

<3-29>  4k-  " , “ii  /6ij.  • 

m x , j=l 

Recall  that  Po  = **0(i)j*l  * So 


’«  ' tl 

r 

a . . « 7 a 6 .... 

ij  s=l  5 °s(l)j 


Thus 


the  multinomial  expansion  of  a 


This  implies  at  once  that  the  expansion  of  (3.29)  contains  a term  of  the  form 


contains  only  positive  summands  from  (2.4)  and  (3.33) 


we  obtain  the  inequality  (3.27) 


Consider  the  identity 


According  to  Lemma  3.1 


From  (3.27)  and  (3.35)  it  follows 


Applying  the  inequality  (3.20)  and  noting  that 


1 we  finally  deduce 


Using  Stirling's  formula  and  the  fact  that  lim  e = 0 we  get 


Hall  (unpublished) . The  problem  is  to  estimate  from  below  the  permanent  of  0-1  matrix 


having  exactly  three  1 in  each  row  and  column.  In  what  follows  we  consider  a larger  class 
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